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Abstract 

The uncertainties in the perturbative and higher twist corrections to the sum rules for Tp^n 
are analyzed. The theoretical predictions for Tp „ are compared with the whole set of ex- 
perimental data and the restrictions on theoretical uncertainties was obtained. Suggestions 
for future experiments are given. 



1 The sum rules for Fp^. Theoretical Status. 

I start with the consideration of the sum rules for the first moments of the spin structure functions 

1 

Tp.nlQ^) = J dxgip,n{x, Q^) (1) 



I will discuss the uncertainties in the theoretical predictions for Fp „ and compare the theoret- 
ical expectations with experimental data. The aim of this consideration is to obtain from the 
experiment the restrictions on the uncertainties in the theoretical description of the problem. 
Consider first the Bjorken sum rule [1]: 

TAQ') - T,m = - - 3.6 [^)^- 20 [^)\ + (2) 

where qa is the axial /3-decay coupling constant and the last term in (2) represents the twist-4 
correction. The perturbation QCD corrections are known up to the third order '^'^^ (there is also 
an estimate of the fourth order term [^]). The coefficients in (2) correspond to the number of 
flavours Nf = 3. 

In what follows I will so often transfer the data to the standard reference point = Ql — 
lO.bGeV'^ - the mean value of Q^, at which EMC and SMC experiments were done. 

Let us discuss the perturbative corrections in (2). Today there is a serious discrepancy in the 
values of ag found from different experiments. The average value of as, obtained at LEP is [^] 

Q;,(m|) =0.124 ±0.007 (3) 
In two-loop approximation this value corresponds to the QCD parameter A3 for three flavours 

A3 = 430 ± lOOMeV (4) 

The data on the T — > hadrons decay give (the first ag correction is accounted) 

Q;^(m^) = 0.178 ± 0.010 (5) 

from which it follows 

A3 = 170 ± 30Mey (6) 

The small error in (6) is caused by the fact that the partial width r(T 3g), from which (5) 
was determined, is proportional to and the as correction to it is small. A strong contradiction 
of (4) and (6) is evident. 

The overall fit '''1 of the data of deep inelastic lepton-nucleon scattering gives in the NLO 
approximation A3 = 250MeV (the error is not given). New data in the domain m| indicate 
lower a,(m|) (SLD I^l : 0.118 ± 0.013,0.112 ± 0.004; OPAL l^l; 0.113 ± 0.012), but the data of 
AMY [101 on e+e- annihilation at = 57.3GeV results in a,(m|) = 0.120 ± 0.005. Finally, 
from T- decay it was obtained ^^^^ : 

as{ml) = 0.33 ± 0.03 (7) 

corresponding to 

A3 = 380 ± 6QMeV 
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But the determination of ag from r-decay can be criticized '^^1 on the grounds that at such a low 
scale exponential terms in (f may persist in the domain of positive — besides the standard 
power-like terms in accounted in the calculation. 

In such a confusing situation I will consider two options - of small and large perturbative 
corrections. In the first case I will take A3 = 200Mel^. Then 

as{Ql) = 0.180 ± 0.010 V^) - ^n{Ql) = 0.194 (8) 

In the second one A3 = AQQMeV and 

as{Ql) = 0.242 ± 0.025 Tp{Ql) - -VM) = 0-187 (9) 

The twist-4 contribution was disregarded in (8) and (9). There is also a discrepancy in its 
value. The value of 6p_„ was determined in the QCD sum rule approach by Balitsky, Braun and 
Koleshichenko (BBK) I^^]; 

6p_„ = -O.OlSGel/^ (10) 

On the other hand, the model t^^l based on connection I^^l of the Bjorken sum rule at large 
with the Gerasimov, Drell, Hearn sum rule at = l^^l gives 

= -O.lSGel/^ (11) 

In what follows I will consider (10) and (11) as two options which correspond to small (S) 
and large (L) twist-4 corrections. I will discuss both approaches of determination of higher 
twist corrections in more details below. From my point of view, no one of these approaches is 
completely reliable and I will use them in comparison with experiments only as reference points. 
I turn now to the sum rules for Fp and F„: 



rp,n(g^) = ^ { 



1_^_3.6(^)2_20(— )3 

TT TT TT . 



{±ga + ^a8)+ (12) 



1-^-0.55(^)2 

OTT TT 



According to the current algebra qa, as and S are determined by the proton (or neutron) 
matrix elements of flavour octet and singlet axial currents 

- 2ms^as =< p, s I JsJ | p, s > - 2ms^E =< p, s | j^J | P, s >, (13) 
-2mSfj,gA =< P, s I jJJ \p,s> (13') 



where 



= 147^7514 + (^7^75^ + 57^755 
In the parton model qa, as and E are equal to 

gA = Au- Ad a8 = Au + Ad- 2AS T, = Au + Ad + As, (14) 

where 

1 

Aq ^ J [q+{x) - q^{x)] q^u,d,s (15) 
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and q± are the quark distributions with the spin parallel (antiparallel) to the proton spin, which 
is supposed to be longitudinal (along the beam). Ag in (12) has the similar meaning, but for 
gluons, as Ag for quarks 
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Ag(Q') ^ jdx [g+{x, Q') - g_(x, Q')] (16) 



Unlike g^, a-s and S which in the approximation used above arc independent and have zero 
anomalous dimensions, Ag anomalous dimension is equal to -1. This means that 

Ag{Q')Q2^^ ~ dng^ (17) 

The conservation of the projection of the angular momentum can be written as 

+ Ag{Q') + L^iQ') = ^ (18) 

where has the meaning of the orbital momentum of quarks and gluons. As follows from 
(17), (18) at high L^iQ'^) must compensate Ag{Q'^), L^iQ"^) ~ —clnQ"^. This means that the 
quark model, where all quarks are in ^-state, failed with increasing. 

Now about the numerical values of the constants gA-, os and E entering eq.(12). gA is known 
with a very good accuracy, gA — 1.257 ± 0.003 [^1. 

Under assumption of SU(3) flavour symmetry in baryon decays is equal to 

as = 2,F-D = 0.59 ± 0.02 (19) 

where F and D are axial coupling constants of baryon /9-decays in SU(3) symmetry and the 
numerical value in the r.h.s. of (19) follows from the best fit to the data t"^^'. The combination 
3F — D can be found from any pair of baryon /3-decays. The comparison of the values of 
3F — D, obtained in this way shows, that the spread is rather narrow I^^l, | Sa^ |< 0.05 and 
at least ag > 0.50. This may be an argument in the favour that SU(3) violation is not large 
here.(See,however, '^^1). It should be mentioned that at fixed Vp^n the uncertainty in ag only 
slightly influences the most interesting quantities E and As. As follows from (12) and (14) 

5E = -5{As) = i(5ag (20) 

ag was also determined by the QCD sum rule method In this approach no SU(3) flavour 
symmetry was assumed and the result 

ag^ 0.5 ±0.2 (21) 

is in agreement with (19), although the error is large. 

What can be said theoretically about E? In their famous paper Ellis and Jaffe assumed 
that the strange sea in the nucleon is nonpolarized. As = 0. Then 

E fa ag fa 0.60 (22) 

This number is in a contradiction with the experimental data pointing to smaller values of E. On 
the other hand, Brodsky, Ellis and Karliner t^^l had demonstrated that in the Skyrme model at 
large number of colours Nc, E ~ 1/A^c and is small. From my point of view this argument is not 
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very convincing: the Skyrme model may be a good model for description of nucleon periphery, 
but not for the internal part of the nucleon determining the value of E (see also P^l). 

An attempt to calculate S in the QCD sum rule approach in the way similar to the ag 
determination ^"^^^ failed. It was found '^^1, that the operator product expansion (OPE) breaks 
down in case of the polarization operator (and/or NNA vertex) for singlet axial current at the 
scale ~ IGeV. (The anomaly was properly accounted in the calculation). This is indicative 
of that the higher dimension operators (instantons?) are of importance in this problem. The 
physical consequence is that one may expect a violation of the Okubo-Zwcig-Iizuka rule in the 
nonet of axial mesons. The similar trouble, perhaps, faces attempts to determine S using the 
so called "U{1) Goldberger-Treiman relation" (for a review see I^^'). So, at this stage, the only 
way to find E is from an experiment, exploiting eq.(12). 

I dwell now on the determination of higher twist corrections. The twist-4 correction was 
calculated by BBK '^^1 , using the QCD sum rule method for the vertex function in the external 
field [26,27] rpj^g result for was given in (10). For bp^n it was obtained 



This result, however, cannot be considered as reliable for the following reasons: 

1. BBK use the same hypothesis as EUis and Jaffe did, i.e. assume that s-quarks do not 
contribute to the spin structure functions and instead of singlet (in flavour) operator consider 
the octet one. 

2. When determining the induced by external field vacuum condensates, which are very 
important in such calculations (sec t^^l), they saturate the corresponding sum rule by ^^-mcson 
contribution, what is wrong. (Even the saturation by r/'-meson would not be correct, since rj' is 
not a Goldstone) . 

3. One may expect that in the same way as in the calculation of S by the QCD sum rule, in 
this problem the OPE scries diverge at the scale ~ IGeV, where the BBK calculation proceeds. 

Even in the case of the Bjorken sum rule, where the mentioned above problems arc absent, 
the value (10) is questionable. was calculated by BBK using the OPE for the vertex 
function and accounting few terms in OPE with operator dimensions from to 8. But the final 
result comes almost entirely from the last term of OPE - the operator of dimension 8. (There are 
other drawbacks in these calculations - see '^^1). It is clear that such situation is unsatisfactory. 
Recently, Oganesian '^^1 had calculated the next term in OPE for 6p_„ determination - the 
dimension 10 term and found that by absolute value it is equal to (10), but has opposite sign, 
so the total result is zero. Of course, we cannot believe in this statement either: it means only 
that the results of the calculations are unstable and the value (10) characterizes the answer by 
the order of magnitude only. 

The other way to determine higher twist corrections suggested in [-"^^1 is based on the con- 
nection of Tp^niQ^) at high with the Gerasimov-Drell-Hearn '^^1 (GDH) sum rule at = 0. 
The idea is the following. For the real photon there is only the spin-dependent photon-nucleon 
scattering amplitude Si{i'), for which we can write an unsubtracted dispersion relation 



Consider the limit i/ — > in (24). According to the F.Low theorem the terms proportional to 
and in the expansion in powers of u of the photon-nucleon scattering amplitude are expressed 
via static characteristics of nucleon. 



6„ , „ = -0.022GeV^ 



(23) 




(24) 
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The calculation gives 

Si{iy),^o = (25) 

where k is the nucleon anomalous magnetic moment: Kp— 1.79, Kn — —1.91. 
Prom (24), (25) the GDH sum rule follows: 

J^G,{u,0)^-l^ (26) 



Till now only indirect check of (26)was performed, where in the l.h.s the parameters of resonances, 
obtained from the nN scattering phase analysis, where substituted. In this way with resonances 
up to 1.8 GeV it was obtained '^^l (see also the second reference in '^"^l and '^^^^) 

l.h.s of (26) r.h.s. of (26) 

proton - 1.03 - 0.803 (27) 

neutron - 0.83 - 0.913 

The l.h.s. and the r.h.s. of (26) are not in a good agreement - a nonresonant contribution is 
needed. The direct check of the GDH sum rule would be very desirable! 

An important remark: the forward spin dependent photon-nucleon scattering amplitude has 
no nucleon pole. This means that there is no nucleon contribution in the l.h.s. of GDH sum rule 
- all contributions come from excited states: the GDH sum rule is very nontrivial. 

In order to connect the GDH sum rule with Tp^^iQ^) consider the integrals [^^^ 

roo fill 

IpAQ')= -G,,p,n{^,Q') (28) 

It is easy to see that at large Q'^ 

IpAQ') = -qT ^pAQ') (29) 

and at (5^ = (28) reduces to the GDH sum rule. For proton Ip is negative at = and 
positive at large Q^, what indicates to large nonperturbative corrections. In l^^l the VDM based 
interpolation model was suggested, describing /p,n(Q^) in the whole domain of Q^. The model 
was improved in the second Refs. 14, where the contributions of baryonic resonances up to 
W = 1.8GeV, taken from experiment, where accounted. The model has no free parameters, 
besides the vector meson mass, for which the value fXy = O.dGeV'^ was chosen. Using this model 
it is possible to calculate the higher twist contributions in (2), (12). The results are presented in 
Table 1, as the ratio of asymptotic F"* with power corrections excluded to the experimentally 
measurable F at given l^^l. 

Table 1. 

Higher twist corrections in GDH sum rule + VDM inspired model. 





2 


3 


5 


10 




1.44 


1.29 


1.18 


1.08 




1.30 


1.20 


1.13 


1.06 


■pas /-p 


1.45 


1.29 


1.18 


1.08 


pas /p 

^ p+nl ^ P+n 


1.47 


1.31 


1.19 


1.08 
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The power corrections, given in Table 1 are essentially larger (except for the case of neutron), 
than the values (10), (23) found in l^^l. It must be mentioned that the accuracy of the model in 
the domain of intermediate Q^, where it is exploited, is not completely certain. 



2 Comparison with experiment. 



When comparing the sum rules (2), (12) with experiment I consider two limiting variants of 
perturbative corrections: small with A3 = 200MeV and large with A3 = 400MeV. (ag is computed 
in 2- loop approximation, it is assumed that the number of flavours Nf = 3). For higher twist 
correction I also consider two limiting options: small (S), given by (10), (23) and large (L), 
determined by the data of Table 1. The contribution of gluons Ag^Q"^) in (12) will be found in 
the following way. Let us assume, that at 1 GeV the quark model is valid and L2(lGeV^) = 
in eq.(18). Taking E=0.3, what is a reasonable average of the data, we have from (18) 



Ag{lGeV^) = 0.35 
The dependence of Ag can be found from the evolution equation 



(30) 



[32] 



Ag{Q') = 



1 + 



2N^ 
bn 



as{Q'^) - a-s 



(/.^)]}a^(/.^) 



+ 



a,. (^2 



where 6 = 11 - {2/2>)Nf = 9, //^ = iGeV^ and E(lG'el/2) 0.3. As the calculation shows, the 
change of scale at which the quark model is assumed to work (say O.SGeV^^ instead of IGeV'^) 
or the of use slightly different S(/x^) in (31) only weakly influence the results for S and As, 
obtained from experimental data. 

I consider the following experimental data (Table 2). 



Table 2. 

The experimental data on Fp . 



Experimental 
group 


The target 


T^p, n) 


Mean Q^{GeV^) 


EMC [33] 


P 


Tp = 0.126 ±0.010 ±0.015 


10.7 


SMC [34] 


P 


Tp = 0.136 ±0.011 ±0.011 


10.5 


E143 [35] 


P 


Tp = 0.127 ±0.004 ±0.010 


3 


E142 [36] 


He' 


r„ = -0.022 ±0.011 


2 


SMC [37] 


d 


Vd = 0.034 ±0.009 ±0.06 


10 






Tp ± r„ = 0.073 ± 0.022 


10 


E143 [38] 


d 


Fd = 0.042 ± 0.003 ± 0.004 


3 






Tp ± r„ = 0.0908 ± 0.006 ± 0.008 


3 
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In comparison with experiment the perturbative and higher twist corrections, as well as 
l^g{Q'^) contributions are calculated for Experimentally, are different in different x-bins 
(higher at larger x). This effect is not accounted in the calculation. The ratio of term to 
a1 term in perturbative corrections is of order of 1 at A3 = AQ{)MeV and Q"^ 2 — SGeV"^ (as 
well as a^/a^ estimate). For this reason we introduce in these cases an additional error equal 
to the . The values of E and As calculated from comparison of experimental data with eq.l2 
are shown in Table 3. (The errors are summed in quadrature). 

Table 3. 

Determination of E and As from experimental data 



Experiment 


A3 


High 


E 


As 




MeV 


twist 








200 


S 


0.21 ± 0.17 


-0.13 ± 0.06 


EMC 


400 


S 


0.29 ± 0.17 


-0.10 ± 0.06 


P 


200 


L 


0.285 ± 0.17 


-0.10 ± 0.06 




400 


L 


0.37 ± 0.17 


-0.07 ± 0.06 




200 


S 


0.30 ± 0.14 


-0.10 ± 0.05 


SMC 


400 


S 


0.39 ± 0.14 


-0.07 ± 0.05 


P 


200 


L 


0.39 ± 0.14 


-0.07 ± 0.05 




400 


L 


0.47 ± 0.14 


-0.04 ± 0.05 




200 


S 


0.28 ± 0.10 


-0.10 ± 0.03 


E143 


400 


S 


0.42 ±0.10 


-0.06 ±0.03 


P 


200 


L 


0.57 ±0.10 


-0.006 ±0.03 




400 


L 


0.71 ±0.10 


0.04 ± 0.03 


E142 


200 


S 


0.60 ±0.12 


0.003 ± 0.04 


n 


400 


s 


0.57 ±0.12 


-0.005 ± 0.04 




200 


L 


0.64 ±0.12 


0.016 ±0.04 




400 


L 


0.61 ±0.12 


0.008 ±0.04 




200 


S 


0.27 ±0.10 


-0.11 ±0.03 


SMC 


400 


S 


0.33 ±0.10 


-0.09 ±0.03 


d 


200 


L 


0.29 ±0.10 


-0.10 ±0.03 




400 


L 


0.34 ±0.10 


-0.08 ±0.03 




200 


S 


0.37 ±0.06 


-0.07 ±0.02 


E143 


400 


s 


0.44 ±0.06 


-0.05 ±0.02 


d 


200 


L 


0.18 ±0.06 


-0.01 ±0.02 




400 


L 


0.54 ±0.06 


-0.015 ±0.02 



Remark: the contribution to E of the term proportional to Agf is approximately equal to 0.06 in 
the case of A3 = 2mMeV and 0.11 in the case of A3 = 400Me\/. 

If we assume that all the analysed above experiments are correct in the limits of their quoted 
errors (or, may be, 1.5 st. deviations), then requiring for the results for E and As from various 
experiments to be consistent, we may reject some theoretical possibilities. A look at the Table 3 
shows that the variant A3 = 400Me\^, L (a contradiction of E143, p and SMC, d results for E) 
and, less certain, the variant A3 = 200M eV, S (a contradiction of E142, n and SMC, d) may be 
rejected. 
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Consider now the Bjorken sum rule. For comparison with theory I choose combinations of the 
SMC data - on proton and deuteron, the E143 data - on proton and deuteron and the E143 data 
on proton and the E142 on neutron i^He). The results of the comparison of the experimental 
data with the theory are given in Table 4. 

Table 4. 

Comparison of the experimental data with the Bjorken sum rule 



Combination 

of exj)eiiiiieiits 


(Tp n) exper. 




High 

twist 








200 


S 


0.193 


SMC, p 


0.199 ±0.038 


400 


s 


0.186 


SMC, d 




200 


L 


0.180 






400 


L 


0.173 


E143,p 




200 


S 


0.182 




0.163 ±0.010 ±0.016 


400 


s 


0.168 


E143,d 




200 


L 


0.145 






400 


L 


0.134 


E143,p 




200 


S 


0.182 


E142 VL{^He) 


0.147 ±0.015 


400 


s 


0.168 


recalculated to 




200 


L 


0.145 


= ^GeV^ 




400 


L 


0.134 



From Table 4 we see again some indications for rejection of variants A3 = 400MeK, L and, less 
certain, A3 = 200Mey, S.(ln the first case there is a contradiction of the theory with the E143,p 
and d data, in the second - with the E143,p , E142, n data). 

At existing experimental accuracy it is impossible to choose from the data the true values of A3 
and twist-4 correction . My personal preference is to the variant A3 = 200Me\^ and to the value of 
twist-4 corrections 3 times smaller than given by the GDH sum rule + VDM inspired model and, 
correspondingly, hp = 0.04 in (12), i.e., 2.2 times larger than the BBK result. The argument in the 
favour of such choice is that at larger A3 there will arise many contradictions with the description 
of hadronic properties in the framework of the QCD sum rules. The recent preliminary SLAC 
data [^^1 on gi{x, Q"^) Q^-dependence indicate that Fp^/Fp = 1 + Cp/Q"^, Cp = 0.25 ± 0.15 what is 
compatible with the estimate above. In this case all experimental data except for E142,n, are in 
a good agreement with one another and the values of E and As averaged over all experiments, 
except for E142,n are 

E = 0.35 ±0.05 As = -0.08 ± 0.02 (32) 

(see also ^''^^ where the values close to (32) were obtained). 

The values of E and As obtained from the E142,n experiment at such a choice of A3 and 
twist-4 corrections, are different: 

E = 0.61 ±0.12 As = 0.01 ±0.04 (33) 

It is impossible to compete (32), (33) by any choice of A3 and higher twist correction. Perhaps, 
this difference is caused by inaccounted systematic errors in the E142 experiment. 
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A remarkable feature of the result (32) (as well as of the data in Table 3) is the large value 
of I As I - the part of the proton spin projection carried by strange quarks. This value may be 
compared with the part of the proton momentum carried by strange quarks 

1// = J dxx [s+{x) + s_{x)] = 0.026 ± 0.006 0.040 ± 0.005 l^^l (34) 

The much larger value of | As | in comparison with contradicts the standard parametrization 

s+{x) + s_(x) = A (1 - xf 

s+{x) - s_{x) = B x-'' {1 - xf (35) 

and requirement of positiveness of s+ and s_, if o; ~ l(pomeron intercept) and 7 < (oi inter- 
cept) . Large | As | and small V2 means that the transitions ss ^ uu + dd are allowed in the 
case of the operator j^s and are suppressed in the case of the quark energy-momentum tensor 
operator 0^j,, corresponding to the matrix element V2. Such situation can be due to nonpertur- 
bative effects and the instanton mechanism for its explanation was suggested Improvement 
of experimental accuracy is necessary in order to be sure that the inequality | As | S> V2 indeed 
takes place. 



3 The calculations of the polarized structure functions 
by the QCD sum rule method 

I mention here only some basing points of the calculations referring for details to [^^1. The 
calculation was performed on the basis of OPE with the account of the unit operator and of 
the square of quark condensate ~ ctg < | ipip \ >^. It was shown that for the bare quark 
loop the Bjorken and Burkhardt-Cottingham sum rules are fulfilled. It was also found that 
for the function gi{x) the results are reliable in a rather narrow domain of intermediate x : 
0.5 < a; < 0.7. In this domain the contribution of li-quarks is much larger than d-quarks, 
gi ^ gf- Therefore, gi ^ (4/9) g^. g^ was calculated and for the mean value of gi in this 
interval it was obtained (at ~ 5 — 10Ge\^^): 

^1(0.5 <x< 0.7) = 0.05 ± 50% (36) 

The large uncertainty in (36) results from the large contribution of nonleading term in OPE and 
from large background at the phenomenological side of the QCD sum rule. The E143 proton '^^1 
and deuteron [^^^ data (the latter under assumption that g"^ is small in this interval of x) give 
roughly the same values: 

^1 (0.5 < X < 0.7) = 0.08 ±0.02 (37) 
This value is in a good agreement with the SMC result t"^^' 

gi (0.4 <x< 0.7) = 0.08 ± 0.02 ± 0.01 (38) 

and with our theoretical expectation (36). 

For the case of the structure function (72 only can be calculated at 0.5 < a; < 0.8, the 
calculation of fails because of large contribution of nonleading terms in OPE. If we assume 
that like in the case of gfi, | 5'2 l"^l fl'2 I) 

^2(0.5 <x< 0.8) = -0.05 ± 50% (39) 
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The E143 data l^^l in this interval of X are: 

52(0.5 <x < 0.8) = 
in a good agreement with (39). 



-0.037 ±0.020 ±0.003 



(40) 



4 What would be desirable to do experimentally in the 
near future? 

1. To increase the accuracy by 2-3 times. 

2. To study the Q^-dependencc (in separate bins in x). 

3. To go to higher (at HERA). The experiment at higher will be informative if only its 

accuracy will not be worse than the existing ones. 

4. To have better data in the domain of small x. (Probably, this also can be done only at 
HERA). 

5. Measurements of two-jets events in polarized deep- inelastic scattering, which will give infor- 
mation about Ag. 

6. To have more data on s-quark distribution in nonpolarized nucleon, particularly at x < 0.1. 

7. To have better data for g2. 

8. To perform a direct check of the GDH sum rule. 
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